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Abstract 

In this paper wc study an approximation scheme for an Hamilton- Jacobi equa- 
tion of Eikonal type defined on a network. We introduce an appropriate notion 
of viscosity solution for this class of equations (see [12]) and we prove that an 
approximation scheme of semi-Lagrangian type converges to the unique solution 
of the problem. 
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1. Introduction 

There is an increasing interest in the study of linear and nonlinear PDEs 

defined on networks since they naturally arise in several applications (internet, 
vehicular traffic, social networks, email exchange, disease transmission, etc.) 
While a theory of linear PDEs on networks is fairly complete (see [10], [11]), 
the study of nonlinear problems is very recent ([6]) and, concerning Hamilton- 
Jacobi equations and control problems on networks, is still at the beginning (see 
[1], [7], [12]). 

It is well known that Hamilton- Jacobi equations in general do not admit reg- 
ular solutions and the correct notion of weak solution is the viscosity solution 
one. Hence all the three papers concerning Hamilton-Jacobi equations aim to 



Email addresses: camilliSdmmm.uiiiroiiial . it (Fabio Camilli), festa@mat.uniromal.lt 
(Adriano Festa), DirkOschieborn.de (Dirk Schieborn) 



Preprint submitted to Elsevier 



December I4, 2012 



extend the concept of viscosity solution to the case of network and, in particu- 
lar, to find the correct transition condition at the internal vertices. But, since 
the papers are motivated by different model problems and therefore they make 
different assumptions on the Hamiltonian at the vertices, the resulting defini- 
tions of viscosity solution are quite different, even if all of them give existence 
and uniqueness of the solution. 

The definition of viscosity solution introduced in [12] satisfies a stability 
property with respect to the uniform convergence. In this paper, we take ad- 
vantage of this property to prove the convergence of a numerical scheme for 
Hamilton- Jacobi equations on a network. For sake of simplicity wc consider 
an Hamiltonian of Eikonal type, i.e. H{x,p) = \p\ — f{x), with a Dirichlet 
boundary condition, but the results can be extended to a more general class of 
Hamiltonians and also to other boimdary conditions. 

Following [4], we introduce a scheme of semi-Lagrangian type by discretiz- 
ing with respect to the time the representation formula for the solution of the 
Dirichlet problem. Wc prove the well poscd-ness of the discrete problem intro- 
ducing an appropriate discrete transition condition and the convergence of the 
scheme to the solution of the continuous problem. It is worth noticing that the 
proof can be adapted to prove convergence of other approximation schemes, for 
example based on finite difference approximation. 

In the second part of the paper we study a fully discrete scheme which gives 
a finite-dimensional problem. The scheme is obtained via a finite element dis- 
cretization of semi-discrete problem. Also for this step of the discretization 
procedure we prove the well posed-ness of the discrete problem and the con- 
vergence of the scheme to the unique solution of the continuous problem. It 
is important to observe that the scheme not only computes the solution of the 
Eikonal equations, but it also produces an approximation of the shortest paths 
to the boundary. 

We also discuss some issues concerning the implementation of the algorithm 
and we present some numerical examples. 

2. Assumptions and preliminary results 

We give the definition of graph suitable for our problem. We will also use 
the equivalent terminology of topological network (see [9] ) . 

Definition 2.1. Let V = {vi, i Cz 1} he a finite colleetion of different points 
in M.-^ and let {tTj, j ^ J} be a finite collection of differentiable, non self- 
intersecting curves in given by 

TTj : [0,lj] M^, Ij > 0, j e J. 

Set Cj := Trj{{0,lj)), Cj := 7rj([0, 1^]), and E := {e^ : j G J}. Furthermore 
assume that 

i) T^j{0),''^j{lj) € V for all j G J, 
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ii) #{ej r\V) = 2 for all j e J, 

in) Bj (lek CV, and #{ej fl e^) < 1 for all j, k G J, j ^ k. 

iv) For all v,w €V there is a path with end-points v and w (i.e. a sequence 
of edges {ej}jLi such that #{ej n e^+i) = 1 and v & ei, w & e^ ). 

Then f := (J^g j C M.^ is called a (finite) topological network in . 

For i G I we set Inci := {j G J : ej is incident tofi}. Given a nonempty set 
Ib C /, we define i9r := {vi, i G Ib} (we always assume i & Ib whenever 
#{Inci) = 1 for some i e /.) We set It--=I\ Ib and T := f \ dT. 
For any function « : F — )• K and each j G J we denote by the restriction of u 
to Cj, i.e. 

:= uo TTj : [0, Ij] M. 

We say that u is continuous in F and write u G C{r) if u is continuous with 
respect to the subspace topology of f . This means that e C{[0,lj]) for any 
j G J and 

u^ttJ'^{vi)) = vl' {'K~'^ (vi)) for any i G I, j,k G InCi. 
We define differentiation along an edge ej by 

d 

dju{x) := djU^{'Kj^{x)) = Q^u\TrJ^{x)), for x G ej, 
and at a vertex Vi by 

d 

dju{x) := djU^^nJ^^x)) = -g^u^ {''^j'^ i^)) for x = Vi, j G InCi. 

Observe^ that tlic; paranictrization of the arcs Cj induces an orientation on the 
edges, which can be expressed by the signed incidence matrix A = {a,j}ij£j 
with 

{1 \i Vi G Ej and 7rj(0) = fj, 
— 1 ii Vi G Cj and nj{lj) = Vi, (2-1) 
otherwise. 

Definition 2.2. Let <fi G C(F). 

i) Let x G Cj, j G J. We say that ip is differentiable at x, if ip^ is differen- 
tiable at 'kJ^{x). 

ii) Let X = Vi, i G It, j, k G Inci, j ^ k. We say that ip is {j, k) -differentiable 
at X, if 

aijdj(fj{TT~'^{x)) + aikdkfkiTTk^ix)) = 0, (2.2) 
where (a-y) as in (2.1). 

Remark 2.1. Condition (2.2) demands that the derivatives in the direction of 
the incident edges j and k at the vertex Vi coincide, taking into account the 
orientation of the edges. 
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We consider the eikonal equation 

\du\ - f{x) =0, a; e r. (2.3) 

where / e C°{T), i.e. f{x) = PinJ^x)) for x e ej, p € C°{[0,lj]), and 
f^{'Kj^{vi)) = f^{T^k^{vi)) for any i G I, j,k G Inci. Moreover we assume that 

f{x) > r? > x&T (2.4) 

Definition 2.3. 

A function u £ USC{T) is called a (viscosity) subsolution of (2.3) in T if the 
following holds: 

i) For any x € Cj , j G J , and for any ip G C{T) which is differentiable at x 
and for which u ~ ip attains a local maximum at x, we have 

\dMx)\ - m ■■= \dj^j{nj\x))\ - fH^jHx)) < 0. 

a) For any x = Vi, i G It, and for any which is {j,k) -differentiable at x 
and for which u — ip attains a local maximum at x, we have 

A function u G LSC{T) is called a (viscosity) supersolution of (2.3) in T if the 
following holds: 

i) For any x G ej, j G J, and for any ip G C(r) which is differentiable at x 
and for which u — ip attains a local minimum at x, we have 

\dMx)\- f{x)>Q. 

a) For any x = Vi, i G It, j G InCi, there exists k G InCi, k ^ j, (which 
we will call i-feasible for j at x) such that for any ip G C{T) which is 
{j, k)- differentiable at x and for which u — ip attains a local minimum at 
X, we have 

\dMx)\-f{x)>o. 

A continuous function u G C(r) is called a (viscosity) solution of (2.3) if it is 
both a viscosity subsolution and a viscosity supersolution. 

Remark 2.2. Let i G It and ip G C(r) be {j, k) -differentiable at x = Vi. Then 

\dMx)\ - fix) = \d,ip,{7r-\x))\ - P{n-\x)) 

= I ± djPk{7T^\x))\ - f{7r^\x)) = \dMx)\ - fix), 

hence in the subsolution and supersolution condition at the vertices, it is indif- 
ferent to require the condition for j or for k. 
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Wc give a representation formula for the solution of (2.3) completed with the 
Dirichlet boundary condition 

u{x) = g{x) xedV (2.5) 

We define a distance-like function 5 : F x F — ^ [0, oo) by 

5(a;,t/):=inf|^ fijis))ds: i > 0, 7 € 

where 

i) 7 : [0, i] ^ r is a piecewise differentiable path in the sense that there are 
to '■= < ti < ■ ■ ■ < tn+i ■= t such that for any m = 0, . . . , n, we have 
7([*m, Wi]) C ej^ for some jm & J, ttJ^ 07 & C^{tm,tm+i), and 



\i{s)\ 



^(^L'°7)(«) 



1. 



ii) y is the set of all such paths with 7(0) = x, ^{t) = y. 

If f{x) = 1, then S{x,y) coincides with the path distance d{x,y) on the graph, 
i.e. the distance given by the length of shortest arc in T connecting y to x. The 
following result is in the spirit of the corresponding results in in [3], [5], [8] 
(for the proof, see [12, Proposition 6.1]) 

Theorem 2.1. Let g iT ^M. be a continuous function satisfying 

g{x) — g{y) < S{y,x) for any x, y G dV. (2-6) 

Then the unique viscosity solution of (2.3) -(2.5) is given by 

u{x) := min{.g(y) + S{y, x) : y G ^F}. (2.7) 

Remark 2.3. It is worthwhile to observe that if supersolutions were defined 
similarly to subsolutions, then the supersolution condition could not be satisfied 
by (2.7). Consider the network F = Uf^^Ci C R? , where ei = {0} x [0,1/2], 
62 = {0} X [—1,0], 63 = [0,1] X {0} and the equation \du\ — 1 = with zero 
boundary conditions at the vertices Vi = (0, 1/2), V2 = (0, —1), U3 = (1, 0). Then 
the distance solution, see Theorem 2.1, is given by u{x) = inf{(i(y,a;) : y € 9F} 
where d is the path distance on the network. The restriction of m to 62 U 63 has 
a local minimum at the vertex vq — (0,0). Hence if tp is a constant function, 
u — if has a local minimum at wq and therefore the supersolution condition is 
not satisfied for the couple (62,63). Instead the arc 61 is VQ-feasible; see the 
definition of supersolution, for both the arcs e-i and 63 . 

3. The approximation scheme 

We consider an approximation scheme of semi-Lagrangian type for the prob- 
lem (2.3)-(2.5). 
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3.1. Semi- discretization in time 

Following the approach of [4] we construct an approximation scheme for 
the equation (2.3) by discretizing the representation formula (2.7). We fix a 
discretization step /i > and we define a function : f M by 

Uh{x) = inf{J-;,(7'') + g{y) : 7'' e B^^^, y e dT} (3.1) 

where J'hil'') = E*Lo ^'/(7m)km| and 

i) An admissible trajectory 7'' = {lm}m=i C F is a finite number of points 
7m = ^jm. (fm) € r such that for any m = 0, . . . , M, the arc 7m7m+i ^ ^jm. 
for some jm e J and := \ hz±lz^\ < 1 

ii) y is the set of all such paths with 7q — x, = y. 

Remark 3.1. Given 7'' € B'^ y, we define a continuous path, still denoted by 
7'', in Bx,y by setting "/'^{s) = Trj^(tm + (^~^^) (t^^-^ - tm)) for s G [mh, (m + 
l)h] «/7m7m+i ^3m- Then, recalling formula (2.7) we approximate 

Mh M .mh M 

f{l{s)ms))\ds = J2 / fhi^Qmlds « J2 hfilt)W 



which shows that (3.1) is an approximation of (2.7). In the continuous case it is 
always possible to assume by reparametrization that [7(5) | = 1. In the discrete 
one we consider instead velocities in the interval [—1,1], since otherwise near 
the vertices the discrete dynamics can move only in one direction. 

Let B(r) be the space of the bounded functions on the network. We show that 
the function Uh can be characterized as the unique solution of the semi-discrete 
problem 

Uhix) = S{h,x,Uh) (3.2) 
where the scheme 5 : M+ x f x B{T) M is defined by 

S{h,x,ip)= inf {(p{xhq) + hf{x)\q\} (3.3) 

qe[-l,l]: Xhqeej 



if a; = Trj{t) G ej 



S{h, X, (f) = inf 



inf {ip{xhq) + hf{x)\q\} 



(3.4) 



if X = Vi, i G It 

S{h,x,ip)=g{x) (3.5) 

if a; e aF 

where, for x = i^j{t), we define Xhq := 'Kj{t — hq). 
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Proposition 3.1. Assume that 

g{x) < inf{J-^(7) + g{y) : 7 e B^^, y e dV} for any x e dT. (3.6) 

Then Uh is the unique solution of (3.2). Moreover uu is Lipschitz continuous 
uniformly in h, i.e. 

\uh{xi) - Uh{x2)\ < Cd{xi,X2) for any xi,X2&f (3.7) 

Proof. Let ui, U2 be two bounded solutions of (3.2) and set Wi{x) = 1 — e^"'(^\ 
for i = 1, 2. Then Wi satisfies 

Wi{x) = S{h,x,Wi) (3.8) 

where 

S{h,x,^)= inf {e-''^(^)l«l¥'(xftg) + l-e-''^(^)l«l} 



S{h,x,(p) = inf 

k^lnci 



if a; = TTj {t) e ej 
inf {e-^f^''^\''\^p{xhq) + 1 - e-''-f(^)l«l} 



ii X = Vi, i & It 

S{h,x,ip) = l-e-5(^) 

if X e ar 

where, for x = i^j{t), Xhq ■= TTj{t — hq). In fact, for any q e [—1, 1] such that 
Xhq € Cj, we have 

^■(^x) = 1 — e~"*^^^ < 1 — e~"*^^''«^~''-''^^^l^l = 1 — e~"*^^''«^e~''-''^^''^' 
1 - e-«'(^''<')) e-'*^(^'l9l + 1 - e-''^(^'l«l = e-''/(^)l«lM;i(xftq) + 1 - e-'^^^^'I'l 

and the first equation in (3.8) follows taking the infimum with respect to q. We 
proceed similarly for the other two equations. 
We have that 

sup \S{h, x,wi{x)) — S{h,x,W2{x))\ < ^sup|wi(a;) — W2{x)\ 
r r 

with /3 = e"'"' < 1, see (2.4). Since 5 is a contraction, we conclude that 
for /i > there exists at most one bounded solution of (3.8) and therefore of 

problem (3.2). 

Now we show the function Uh is a bounded solution of (3.3)-(3.5). It is 
always possible to assume, by adding a constant, that g > 0. It follows that 
Uh>0. Moreover it is easy to see that 

Uh{x) < ||/||ooSupd(a;,ar) + sup g{x). 

xer xedr 
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To show (3.5), observe that we have Uh{x) ^ g{x) for x G dV if and only if there 
is some z G dT such that g{x) > g{z) + J^hil'^) for some 7'' S B'^ ^ which gives 
a contradiction to (3.6). 

We consider (3.3) and we first show the "<"-inequahty. For x € ej and for 
q e [—1, 1] such that Xhq € Cj, let y & dT and 7^* € B^^ y be e-optimal for 
Uh{xhq). Define 7'* = {7,'*}Lo with 7^ = x, 7^ = x,,,. Hence 7^ U 7'^ e B^^^ 
(with x/jq counted only one time in U 7'') and 

«/.(a;) < g{y)+Fh{l''^li) < 9iy)+J'h{l^) + hf{x)\q\ < Uh{xhq) + e + hf{x)\q\. 
To show the reverse inequality, assume that for some a; G F, 

Uh{x) < ^ inf {uh{xhq) + hf{x)\q\} - 5. 

for 5 > 0. Given e < 5, let y € 9F and 7^' = {7^}^^=o ^ be e-optimal for 
X. By the inequality 

9{y) + ^hhxy) - e < Mft(a;) < Uh{xhq) + /i/(a;)|g| - 5 

it is clear that if y = x/j^ for some q G [—1,1] we get a contradiction. Define 
1^ = lx,v \ 1^ where 7^ = {7f}i=o with 7^ = x, 7^ = a;/,g. Since 7'' := 
7^,j,\7'^eB^^^,^ wehave 

9{y) + ^/^(7'') = 9{y) + J'hiliy) - Fhil^) < Uhixhq) + e-5 

a contradiction to the definition of Uh and therefore (3.3). The equation (3.4) 
can be proved in a similar way. 

We finally show that the function Uh is Fipschitz continuous in F, uniformly in 
h. Consider first the case of two points in the same arc, i.e. x'l, X2 G e.j for 
some j e J. Given e > 0, denote by 7'' = {7^} € ■B^i,x2 




(3.9) 

where \iT~^{jra)-T^~^{lm+i)\ < /i for TO = 0, . . . , M. Let y e ar and 7^ e B'^^ y 
be e-optimal for Xi. Then 7^ U 7^* G B^^ y and 

Uh{x2) < 9{y) + Fhili U 72^) < 5(2/) + J'hili) + J^/.(72 ) 

M 

< Uh{xi) + C' XI 'ikj(*m-i-i - '^j{tm)\ + e < M/»(xi) -|- Cd(xi,X2) -|- 2e 

TO=0 

Exchanging the role of xi and X2 we get 

\uh{xi) - Uh{x2)\ < Cd{xi,X2) (3.10) 
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If xi,X2 G r, let 7 be such that |7(s)|(is < .T2) + e and {ej^}^^i C J 
such that 7([0,T]) C ^m=i^jm- For each one of the couples (xi, t;^^), (fj„,fjv„+i) 
for m = 1, . . . , M and , 0:2) define a trajectory 7^ as in (3.9). Then define 
7'' e B^,,,, by 

a;i, fc = 0; 
X2, m = M. 

where M = E^tJ' M,. For tk = 7r7^^(7fe^), fc = 1 ^i-i, • • • , E™ 1 + 
Mm — 1, then we have tk+i — tk = hqk with \qk\ < 1. Let y G dV and 7^ € -6^1,2/ 
be e-optimal for xi. Then 7{' U 7'' G i?^2,y and 

«?»(X2) < 9{y) + Mil u 72 ) < 9{y) + + ^h{.ii) 

M 

< Uh{xi) + /i|gfc|/(7fe) + e < Uh{xi) + Cd(a;i, 0:2) + 2e. 
fe=o 

Exchanging the role of xi and a;2 we get (3.10) ■ 

Remark 3.2. By Remark 3.1 and the continuity of f, assumption (2.6) implies 
g{x) < inf{jr^(7) + g{y) : 7 e B^^, y G dT} + Ch for any x,y G dT. 

Moreover, if g = on dT, the condition (3.6) is satisfied since J^hil^) > for 
any 7''. 

Theorem 3.1. Assume (3.6) for any h > and (2.6). Then for h 0, the 
solution Uh of (3.2) converges uniformly to the unique solution u of (2.3)-(2.5). 

Proof, we first observe that (2.3) can be written in equivalent form as 

sup {-qdu{x) - f{x)\q\} = 
?e[-i,i] 

By (3.7), Uh converges, up to a subsequence, to a Lipschitz continuous function 
u. Wc show that u is a solution of (2.3) at x E T. We will consider the case 
X = Vi E It, as otherwise the argument is standard (see f.e. [2. Tli.VI.1.1]). 

To show that w is a subsolution, choose any j, k G Inci, j ^ k, along with an 
(j, A;)-test function 99 of u at x. Observe that it is not restrictive to consider x to 
be a strict maximum point for u — tp, since we otherwise consider the auxiliary 
function (ps{y) := (f{y) + 5d{x,yY for 5 > with dm{d{x, [x)) = for 

m = j and m = k. Then there exists r > such that u — (p attains a strict local 
maximum w.r.t. Br{x) at x, where Br{x) {y G F : d{x,y) < r}. Moreover x 
is a strict maximum point for u — <^ also in B := Br{x) fl (cj U Cfe). Now choose 
a sequence w/^ — )• for /i — )• with 

suY)\u{x)-Uh{x)\<uJh (3-11) 
r 
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and let yu be a maximum point for u/j — <^ in B. Up to a subsequence, yh 

z € B. Moreover, 

u{x) - ip{x) -uJh< Uh{x) - ip{x) < Uh{yh) - fiUh) < u{yh) - <p{yh) + f^h- 

For /i — ^ 0, we get u{x) — (p{x) < u{z) — (p{z)- As a; is a strict maximum point, 
we conclude x = z. Invoking 

u{x) + (p{yh) - (p{x) -ojh< Uh{yh) < u{yh) + uih 
we altogether get 

lim yh = X, lim Uhivh) = u{x) (3-12) 
We distinguish two cases: 

Case 1: yh ^ X. Then yh G e„i with cither m = j or m = fc. Since Uh — ^ 
attains a maximum at yh, then for yh — iTmith) and yhq = iTmith — hq) € em 

and therefore 

(_ ^{-;n\ykq))-HnM) _ hr(^y^)lql\ < (3.13) 

<?6[-l,l]:yh,6e™ I " J 

The set {g e M : 7r„(f — hq) G e^} contains for h small enough either [—1, 0] if 
ai,m = 1 or [0, 1] if Oi,,,, = —1. Passing to the limit for /i ^ in (3.13), since 
f^{x)\q\= nx)\-q\weget 

sup {qdmifix) - f{x)\q\} <0. 

96[-l,l] 

Case 2: yu = x. Since Uh — f attains a maximum at x, then for x = nj{th) 
and yhq = nj{th - hq) € ej 

UhiVh) - viVh) > u{yhq) - ViVhq) 

and therefore 

sup J_ ^i(w-^i(./.) _,^.(^^)|^|Uo 

qe[-l,l]:j/^,eej ^ " J 

The set G K : ■JTj{t — hq) G ej} contains for h small enough either [—1,0] if 
ttij = 1 or [0, 1] if aij = —1 and passing to the limit for h ^ we conclude as 
in the previous case that 

sup {qdjip{x) - f{x)\q\} <0. 

To show that u is a supersolution, we assume by contradiction that there exists 
j G Inci such that for any k G Inci, k ^ j, there exists a {j, fc)-test function 
of « at a; for which 

sup {qd.Mx)- f{x)\q\} <0. (3.14) 
9e[-i,i] 
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By adding a quadratic function of the form —a^d^x, y)^ to the function (pk we 
may assume that there exists r > such that u — ipk attains a strict minimum 
in Br{x) at x. Observe that x is a strict minimum point oi u — ifk also in 

Since for any h, there exists such that 

<K)= , inf _ {u1'Mt-hq))+hfHviM} 

we may assume, up to a subsequence, that there exists k G Inci such that kh = k 
for any ft- > 0. 

Let Uh be a minimum point of Ut — fk in Bfc and let LOh be as in (3.11). As in the 
subsolution case, we prove that (3.12) holds. If yh ^ x, we have for yh = nm{th) 
and th — hq £ 

UhiVh) - viVh) < u{T^m{th - hq)) - (piwmith - hq)) 

and therefore 



sup l ^tMt,-H))-^l^{y,) _ ^ ^ ^ 

<?6[-l,l]: -n-^Ct-Meem I 

for either m = j or m = k. li yn = x, wo got 



sup _ ^,Mt,-hq))-^iiy,) _ ^^,^^^1^1 1 ^ Q 

ge[-l,l]:7r,(t-ftg)ee,- [ J 

Arguing as in the subsolution case we get for /i — )• 

sup {qdjip{x) - f{x)\q\} > 0. 

96[-l,l] 

which is a contradiction to (3.14). 

We conclude the proof by observing that the uniqueness of the solution to 
(2.3) implies that any convergent subsequence Uh must converge to the unique 
solution u of (2.3)-(2.5) and therefore the uniform convergence of all the se- 
quence Uh to u. m 

3.2. Fully discretization in space 

In this section we introduce a FEM like discretization of (3.2) yielding a fully 
discrete scheme. For any j G J, given Ax-' > we consider a finite partition 

={t{ = o<...<ti<---<ti,^=ij} 

of the interval [0,lj] such that IP-'] = maxi_,,.^Mj(i^ — t^n-i) — ^x-' . We set 

Ax = max Ax-'', M = ^Mj (3.15) 
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The partition induces a partition of the arc ej given by the points 
xL = '^j{iL), m = l,...,Mj. 

and we set X^x = ^jeJ ^mLi xj^- 

In each interval [0, Ij] we consider a family of basis functions {,94}m=o 
space of continuous, piecewise linear functions in the intervals of the partition 

. Hence are piecewise linear functions satisfying P^{tk) = 5mk for m,k 
{!,..., Mj} < /34(t) < 1, Efli PLit) = 1 and for any t e [0, Ij] at most 2 
are non-zero. We define /3j : ^ R by 

Given W G M*^ we denote by Iax [W] the interpolation operator defined on the 
arc Sj by 



Mi Mi 



X € Cj. 



We consider the approximation scheme 

U = S{Ax,h,U) (3.16) 
where the scheme <S = {<S(Aa;, h, W)}j^j is given by 

Si,{Ax,h,W)= inf {I^[W]{xUq)) + hfixi)\q\} (3.17) 

if e Bj 

Si,{Ax,h,W)= inf {l''[W](xtiq)) + hfixt)\q\} (3.18) 

<!£[-l,ll:x^(g)eefc 

if x^^ = ViG It 

Si,{Ax,h,W)=givi) (3.19) 

if a;^ = Vi,iG Ib 

for xi^{q) = Tr^{ti^-hq). 

Proposition 3.2. For any Ax > with Ax < h/2, there exists a unique 
solution U e to (3.17)-(3.19). Moreover, defined UhAx{x) = Iax[U], if 
Ax = o{h) for h ^ 0, then UhAx converges to the unique solution u of (2.3)- 
(2.5) uniformly in T. 

Proof. We show the boundedness of a solution to (3.16) by induction. For 
this purpose we number the nodes Xi such that (i(a;j+i,9r) > d{xi,dT) for all 
i = 1, . . . , M, and claim that 

\Ui\< sup \g{x)\ + h{Lg + Mf) + 2Mfdixi,dT). 

xedT 
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For each xi with d{xi,dT) < h this estimate is immediate. Now assume the 
assertion is true for all Xi with i = 1, — 1. For xi G ej by (3.16) we obtain 
the inequality 

Ui < hf{xi)\q\ +F[U]{xi{q)) < hMf +F[U]{xi{q)) 

for any g e M" with < 1 and xj (q) G ej. Choosing q such that d{xj {q), dT) = 
d{xi,dT) — h and using Aa; < h/2 we obtain that the value X-' [[/](a;; (g)) only 
depends on nodes Xi^ with d{xi^,dr) < d{xi,dT) — h/2, thus ik < I- Picking that 
node ajjj. such that i7,^ becomes maximal, and using the induction assumption 
we can conclude 

Ui < Mfh + < Mfh + sup \g{x)\ + h{Lg + Mf) + 2Mf{d{xi,dT) - h/2) 

xedr 

i.e. the assertion. 

To show existence of a unique solution U we apply the transformation 

to (3.16). Hence f7 is a solution to 

U = §{Ax,h,U) (3.20) 

where 

S^^{Ax,h,W) = inf {e-''/(^-)j^ [VF](a;^„(g)) + 1 - e-''-^(^-)l«l} 

ge[-l,l]:x^(g)eej 

if xi^ e ej 

Si{Ax,h,W)= inf {e-^-f(^-)l'=[W](a;^(g)) + l-e-''-f(^-)l«l} 

<;e[-i,i]: x^(9)eej. 

if xl^ = Vi€ It 

Sl^{Ax,h,W) = l-e^s^"') 

if x{^ = Vi, i & Ib 

As in the proof of Proposition 3.6 we show that 5 is a contraction in M*^ and 
we conclude that there exists a unique bounded solution to (3.20) and therefore 
to (3.16). 

To show the convergence of w^ax to u, we set U/j = 1 — e""'*, UhAx = 
1 — e~"'^^^ and we estimate for x e Cj 

\uh{x) - UhAx{x)\ < \uh{x)-P[U'']{x)\ + \P[U'']{x)-P[U]{x)\ (3.21) 

where U arc the vectors of the values of Uh, UhAx at the nodes of the grid. 
By the Lipschitz continuity and boundedness of Uh we get 



\uh{x)-X^[U^]{x)\ < CAx (3.22) 
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with C independent of h. Moreover, by (3.8) and (3.20) we get for Xk = 
ir~^{tk) e Cj, Xhq ■■= TTj{tk - hq) and since x'^{q) = x^q 

|C/fc^ - tJk\ < e-''^^-''^\uh{xhq)-V[U\{xi{q))\ < e-'"'||% - %ax||oo (3.23) 
where r] as in (2.4). Substituting (3.22) and (3.23) in (3.21) we get 

C 

and therefore, taking into account Theorem 3.1, we have that if Aa; = o{h) for 
/i — ^ 0, then UhAx converges to u uniformly on F. ■ 

4. Implementation of the scheme and numerical tests 

In this section we discuss the numerical implementation of the scheme de- 
scribed in the previous section and we present some numerical examples. We 
remark again that the most interesting feature of our approach is that it is 
intrinsicaUy one-dimensional, even if the graph is embedded in M.^ . For this 
reason it does not present the tj^pical curse of dimensionality issue which is 
usiiaUy encountered in solving Hamilton- Jacobi equations on M.^ . 

The numerical implementation of semi-Lagrangian schemes has been exten- 
sively discussed in previous works (see for example the Appendix B in [2]), hence 
the only regard is due to vertices, where the information could come from dif- 
ferent arcs. We briefly describe the logical structure of the algorithm we use to 
compute the solution. 

Let A be the m, x m incidence matrix defined in (2.1). We also define a ma- 
trix BC which contains the information on boundary vertices, in particular: 
BC{-, 1) represents a boundary vertex and BC{-, 2) = the value of the Dirichlet 
datum at that vertex. The number of the edges is at most n = (™~^)™ and, 
after having ordered the edges, we define the auxiliary edges matrix B G M^'" 
where the i-row contains the following information: 

• B{i,l) = #knot where the i-arc starts, 

• B(i, 2) = :/^knot where the i-arc ends, 

• B{i,3) = length of the discretized i-arc. 

We choose the same discretization step Aa; = Axi for every edge, so that the 
approximated length of the edge i is Li = trunc{^^^^) G N+ and we consider 
a finite partition 

P' = {4 = 0, t{ = Ax, 4 = 2Aa;, • • • , t\,^_, = (M^ - l)Ax, 4. = Bit, 3)} . 

(4.1) 

The matrix C, contains the grid points of the graph, i.e. for the edge i 

C{i,j) = ni{ti) j = 0,...,Mi (4.2) 
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Finally, wc denote by U{i,j) the the approximated solution at the point C{i,j) 
point. We solve the problem using the following iteration 



HJ-networks algorithm. 



1. Initialize 

U = Uo; 
it=0; 

2. Until convergence, Do 

3. for i=0 to n 



4. If there is an s s.t. B{i, 1) = BC{s, 1) 

5. then C/(i,0) = BC(s, 2); 

6. else 

7. U{i,0) = mm{mm^k\A{B{i,i),k)=i} {l[U]{C{k, ^))} , 

min{fe|^(B(.,i)=-i} {l[U]{C{k, B{k, 3) - ^))}} + hf{C{i,j)) 

8. for j = to S(i,3) - 1 

9. U{i,j) = min„e[-i,i] {l[U]{C{i,j + ^))} + hfiC{i,j)) 

10. If there is an s s.t. B(i, 2) = BC{s, 2) 

11. then U{i, B{i, 3))=BC(s,2); 

12. else 

13. U{i,B{i,3))=min{mm^k\AiBi^.2)=l}{I[U]{C{k,^))}, 
min{,|^(B(,,2)=-i} {l[U]{C{k,B{k,3) - ^))}} + hf{C{i,j)) 



14. re-initialize vertex on U 

15. EndDo 



The interpolation I[U]{C{i,x)) is the usual linear interpolation, i.e., said 
t{x) = trunc{x) 

I[C]{x) = C{i,t{x)) + ^-^^^^[C{i,t{x) + 1) - C{z,t{x))] 

numi,^)) = unA^)) + am - cii,tix))) ^;^l'^z^^ 

(4.3) 

Remark 4.1. The order given to the edges, which is necessary to define the 
previous iteration, brings some additional problems that we have to consider: 

• At the end of each iteration of the method, the values of the solution at a 
same vertex, which is contained in different arcs, could be different. Hence 
we make a re-initialization, choosing for every vertex the minimum of the 
previous values. 
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• It is also important that the initial guess Uq of the solution we use to 
initialize the algorithm is greater than the solution. In fact, if this con- 
dition is not satisfied, for particular choices of the discretization step the 
algorithm could generate a nan correct minimum. 

Scheme ofihe Network 
1 r ^-^ • 

0.8 ■/ \ / \ 

0.6 - / \ 1 \ 

Q.4 ■ / \ / \ 

D.3 / \ I \ 

Oil \ II \ • 

-0.2 ■ \ \ 

-0.4 - \ / \ / 





Arcs 

• Knots 
■ BoundKnt 





Figure 1: Test 1, structure of the graph. 



In the first test we consider a five knots graph with two straight arcs and 
two sinusoidal ones (see figure 1). The only boundary knot is the one placed at 
the origin and the value of the solution at this knot is fixed to zero. The cost 
function is constant, i.e. f{x) = 1 on F. In this case the correct solution is 

u{x) —dist(x,0) ~ \x2\ for the straight arcs 

/•kil (4.4) 

u{x) — / (■\/l + (2tt cos27rt))dt for sinusoidal arcs 

"'0 

An approximated solution is shown in Figure 2. In Table 1, we compare the 
exact solution with the approximate one, obtained by the scheme. We observe 
a numerical convergence to the correct solution in i2-norm and in the uniform 
one. As uniform norm we consider the maximum of the uniform norm of the 
error on every arc and as i2-iiorm the maximum of the L2-iiorm on every arc. 
We can observe an order of convergence close to 0.5 that is the typical theoretical 
order of convergence in the uniform norm of semi-Lagrangian schemes in R", 
(see for instance [4]). 
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Approjimsted Value Function 




Figure 2: Test 1, Ax = 0.025. 



Ax = h 


ll-lloo 


Ord(Loo) 


II-II2 


Ord{L2) 


0.2 


0.1468 




0.1007 




0.1 


0.0901 


0.7043 


0.0639 


0.6562 


0.05 


0.0630 


0.5162 


0.0491 


0.3801 


0.025 


0.0450 


0.4854 


0.0402 


0.2885 


0.0125 


0.0321 


0.4874 


0.029 


0.4711 



Table 1: Test 1. 



Scheme ofihe Network 




-1 -0,6 -O.B -0.4 -0.2 0,2 0.4 0.6 0.0 1 



Figure 3: Test 2, structure of the graph. 
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Figure 4: Test 2, Ax = 0.1. 



In the second test we present a more complicated graph with two boundary 
vertices and a several connections among the arcs. Also in this case, we consider 
a constant cost function f{x) = 1 on F. In Table 2 and in Figure 4 we show our 
results. 

In this case we observe an improvement of order of convergence with respect to 
the previous example. This is due to the fact that the graph is compose of only 
straight arcs and this reduces the error due to the piecewise linear discretization 
of the arcs. 



Ax^ h 


ll-lloo 


Ord{L^) 


II-II2 


Ord{L2) 


0.2 


0.1716 




0.0820 




0.1 


0.0716 


1.2610 


0.0297 


1.4652 


0.05 


0.0284 


1.3341 


0.0127 


1.2256 


0.025 


0.0126 


1.1611 


0.0072 


0.8188 


0.0125 


0.0056 


1.1699 


0.0037 


0.9605 



Table 2: Test 2. 
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In the third test we consider a five knots graph (figure 5), with a running 
cost which is not constant. For any point on the graph x — (xi,X2) G F, we 
take f{x) = 10(xi — I) + i], hence f{x) > > for x £ T. In the example, 
we set Tj = 10"^". The graph of the approximate solution is shown in Figure 6. 
Also in this case we provide a experimental table of convergence for the error 
(Table 3). In absence of an exact solution we compare the approximation for 
various grid sizes with a discrete solution Uex on a fine grid (Aa; = 0.005). 



Ax ^ h 


ll-lloo 


Ord(Loo) 


II-II2 


Ord{L2) 


0.2 


0.3800 




0.2078 




0.1 


0.1800 


1.078 


0.0855 


1.2812 


0.05 


0.08 


1.1699 


0.0419 


1.029 


0.025 


0.035 


1.1926 


0.0222 


0.9164 


0.0125 


0.0166 


1.0762 


0.0103 


1.1079 



Table 3: Test 3. 




Figure 7: Test 4, structure of the graph. 

As our last test we consider a graph with several boundary points and a 
more complicated running cost function /. A representation of this graph is 
shown in Figure 7. We consider the following function / 

f{xi,X2) = 2.1 - sin(47rxi) + cos(47ra;2) (4.5) 

obviously, because of the regularity of this function, its restriction on the arcs 
of the graph is continuous. 

In Table 4 we show a comparison for the error in various grid steps. Also in this 
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Figure 8: Test 4, Ax = 0.05. 



case, in absence of the correct solution, we consider as correct the approximation 
on a fine grid (Ax = 0.005). 



Ax ^ h 


ll-lloo 




II-II2 


Ord{L2) 


0.2 


0.7049 




0.3676 




0.1 


0.2925 


1.2690 


0.1557 


1.2394 


0.05 


0.1460 


1.0025 


0.0777 


1.0028 


0.025 


0.0728 


1.0040 


0.0320 


1.2798 


0.0125 


0.0375 


0.9570 


0.0108 


1.5670 



Table 4: Test 4. 
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